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We suggest a time-dependent dynamical mean-field-hydrodynamic model for the collapse of a 
trapped boson-fermion condensate and perform numerical simulation based on it to understand 
some aspects of the experiment by Modugno el al. [Science 297, 2240 (2002)] on the collapse of 
the fermionic condensate in the 40 K- 87 Rb mixture. We show that the mean-field model explains 
the formation of a stationary boson-fermion condensate at zero temparature with relative sizes 
compatible with experiment. This model is also found to yield a faithful representation of the 
collapse dynamics in qualitative agreement with experiment. In particular we consider the collapse 
of the fermionic condensate associated with (a) an increase of the number of bosonic atoms as in 
the experiment and (b) an increase of the attractive boson-fermion interaction using a Feshbach 
resonance. Suggestion for experiments of fermionic collapse using a Feshbach resonance is made. 

PACS numbers: 03.75.-b, 03.75.Nt 
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I. INTRODUCTION 

Recent successful observation of condensed boson- 
fermion mixtures of trapped alkali atoms by different 
experimental groups 0, U 0, has initiated the inten- 
sive experimental studies of different novel phenomena 
Among these experiments there have been stud_- 
ies of condensate of two components of 40 K [l| and 6 Li 
atoms. Condensation of boson-fermion mixtures 6,7 Li [_ 
23 Na- 6 Li and 87 Rb- 40 K [10] have also been reported! 
The collapse of fermionic condensate in a boson-fermion 
mixture 87 Rb- 40 K has been observed and studied by 
Modugno et al. p| . In this paper we consider the collapse 
dynamics of fermions in a boson-fermion mixture us- 
ing a coupled time-dependent mean-field-hydrodynamic 
model. The bosonic part is treated by the mean-field 
Gross-Pitaevskii (GP) equation and the fermionic part 
is treated by a hydrodynamic model. The present mean- 
field-hydrodynamic model could be considered to be a 
time-dependent extension of a time-independent model 
suggested for the stationary equilibrium states by Ca- 
puzzi et al. [9j- 

Because of the Pauli exclusion principle the fermions in 
spin parallel states in the condensate experience a strong 
repulsion. This is the dominating interaction at short dis- 
tances and avoids the collapse of a fermionic condensate. 
In contrast an attractive bosonic condensate larger than 
a critical size is not dynamically stable [ljj. However, 
if such a bosonic condensate is "prepared" or somehow 
made to exist it experiences a dramatic collapse and ex- 
plodes emitting atoms leading to a condensate of small 
size with a high central density. Under high pressure 
three-body recombination takes place with the emission 
of energy leading to explosion and loss of atoms. The 
possibility of collapse was first suggested and observed 
in 7 Li atoms [l0|. A dynamical study of controlled col- 
lapse and explosion has been performed by Donley et 
al. on an attractive 85 Rb bosonic condensate, where 



they manipulated the inter-atomic interaction by chang- 
ing the external magnetic field exploiting a nearby Fesh- 
bach resonance In the vicinity of a Feshbach reso- 
nance the atomic scattering length a can be varied over a 
huge range by adjusting an external magnetic field, thus 
transforming a repulsive bosonic condensate to a highly 
attractive one. There have been many theoretical pllll4| 
studies to describe different features of the experiment by 
Donley et al. 

In the classic study of a trapped boson-fermion mix- 
ture Modugno et al. [5j simulated a situation of a strong 
effective interaction between fermions via a strongly at- 
tractive boson-fermion interaction which may play an im- 
portant role in the collapse of the fermionic condensate. 
The presence of a strong attractive interaction between 
bosons and fermions can induce an effective attraction 
between the fermions 0|. If this overall effective at- 
traction among the fermions can overcome the Fermi 
repulsive pressure there is a possibility of the collapse 
of a fermionic condensate in a boson-fermion mixture. 
Modugno et al. [{| showed in an experiment with 87 Rb- 
40 K mixture that it is indeed the case and produced 
some results of the collapsing dynamics of this system 
which we would like to understand using a mean-field- 
hydrodynamic model. 

Previously, in addition to the study of the collapse 
of a p urely boson ic condensate 0, fljl> we a l so investi- 



gated [lg, [171, [18| the collapse dynamics in a two-species 
bosonic condensate initiated by an attraction between in- 
terspecies bosons using the mean-field coupled GP equa- 
tion It was found that even when the intra-species 
interaction is repulsive one can have a collapse in one or 
two species due to the interspecies attraction which sim- 
ulates an effective intra-species attraction. The situation 
is qualitatively similar in the present boson-fermion mix- 
ture. The boson-fermion attraction in a boson-fermion 
mixture can induce an attraction in the fermionic sys- 
tem leading to collapse. However, there is no coupled GP 
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equation in the case of boson-fermion mixture which com- 
plicates a theoretical description. Instead, in this paper 
we use a time-dependent version of a recently suggested 
equations of generalized hydrodynamics . These equa- 
tions allow us to span the entire range of boson-fermion 
interaction as in the coupled GP equation for bosons . 

The equilibrium properties and the phase diagram of 
a boson-fermion mixture have been studied by several 
authors 0] to obtain good agreement with experiment 
using a mean-field-type description [2(j . All of them em- 
ployed a time-independent formulation. To the best of 
our knowledge the study of nonequilibrium properties of 
fermionic collapse using a time-dependent formulation is 
considered for the first time in this paper. Our findings 
arc in agreement with the experiment by Modugno et al. 
i. 

In Sec. II we present our time-dependent mean-field 
model. This consists of a set of coupled partial differen- 
tial equations involving the bosonic and fermionic con- 
densate wave functions. We introduce an appropriate 
loss term due to three-body boson-fermion recombina- 
tion. In Sec. Ill we present our results for stationary 
boson-fermion wave functions as well as fermion collapse 
initiated by a growth in boson number and boson-fermion 
interaction. Our stationary results are consistent with 
experiment p| and other numerical studies 0, 0, |20| . 
The present study also yields a faithful representation 
of the collapse dynamics as observed experimentally [l| . 
Finally, a summary of our findings is given in Sec. IV. 



II. NONLINEAR 
MEAN-FIELD-HYDRODYNAMIC MODEL 

The time-dependent Bose-Einstein condensate wave 
function ^(r;^) at position r and time t allowing for 
atomic loss may be described by the following mean-field 
nonlinear GP equation |8( 
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9 B (r;t) = 0, 

(2.1) 



with normalization J dr|5 , s(r; t)\ 2 = N B . Here mj is the 
mass and Nb the number of bosonic atoms in the con- 
densate, jib = |5 , s(r;t)| 2 is the boson probability den- 
sity, g BB = ^n 2 a BB /m B the strength of inter-atomic 
interaction, with cibb the scattering length. The trap 
potential with spherical symmetry may be written as 
Vb(t) = ^mBU B r 2 where ujb is the angular frequency. 
The probability density of an isolated fermionic conden- 
sate n_F = |\E'F(r)| 2 is given by 
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where ^^(r) is the condensate wave function, A = 
?i 2 (67r 2 ) 2 / 3 /(2mF), £f is the Fermi energy, and mp is 
the fermionic mass. The spherical trap is given by 



Vf(i") = \m,F'jJpT 1 . The number of fermionic atoms Np 
is given by the normalization J dr\^p(r)\ 2 = Np- 

However, it has been suggested that if there is an in- 
teraction between the fermions and bosons, Eqs. 12. 1|) 
and get modified to [H 
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where g B p = 2Trh 2 a B p/mR with the boson-fermion re- 
duced mass mpi = m B mp /(m B + mp) where mp is 
the mass of fermionic atoms. The interaction between 
fermions in spin polarized state is highly suppressed due 
to Pauli exclusion principle and has been neglected in 
Eqs. i|2.3|) and i|2.4|) and will be neglected through- 
out this paper. Also, we shall always assume boson- 
fermion attraction (a B F < 0) and boson-boson repulsion 
iflBB > 0). 

Capuzzi et al. |9| developed a set of time-independent 
equations for boson-fermion condensate from a consider- 
ation of hydrodynamic motion of two condensed fluids, 
one bosonic and another fermionic, in a spherical trap at 
zero temperature. The interaction between bosons and 
between bosons and fermions are described by contact 
potentials parametrized by coupling constants gBB and 
gsF defined above. They derived the following mean- 
field energy functional for the system 
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The first integral on the right-hand side of Eq. (|2.5(l is the 
Gross-Pitaevskii energy functional of the bosons and is 
related to the (nonlinear) Schrodinger equation • How- 
ever, the second integral, although bears a resemblance 
with the first, is derived from the hydrodynamic equa- 
tion of motion of the fermions and is not related to a 
Schrodinger- like equation Hence, the derivative term 
in the second integral has a different mass factor 6771^ and 
not the conventional Schrodinger mass factor 2m B as in 
the first integral. Finally, the last integral corresponds 
to an interaction between bosons and fermions. 

Equation (|2.5() corresponds to the following La- 
grangian density 



2 ^ 

i=B,F 



dt 



dt 



3 



h 2 \^ B \ 2 

2m B 

h 2 \V^ F \ 2 
6m f 



V b \^b\ 2 + ^9bb\^b\ 4 



V f \^ f \ 2 + -A\^ f \ 10 / 3 
5 



9bf\^f\ 2 \^b\ 2 - 



(2.6) 



The mean-field dynamical equations for the system are 
just the usual Euler-Lagrange equations 
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where Xk,k — 1,2,3 are the three space components, 
i = B corresponds to the bosonic wave function and 
i = F corresponds to the fermionic wave function. With 
the Lagrangian density l|2.6[) these equations of motion 
become: 
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The normalization of the wave functions are given by 
f dr\^ B (r;t)\ 2 = N B and / dr|^ F (r;<)| 2 = N F . It is 
understood that Eq. (|2.4|l results as a consequence of 
Thomas- Fermi- type approximation to Eq. I|2.9|l . when 
the kinetic energy term in the latter (equation) is ne- 
glected at low temperatures. 

Equations l|2.8|l and l|2.9Jl are the desired time- 
dependent mean-field-hydrodynamic equations. We shall 
use these equation for the study of the collapse of the 
fermionic atoms in a coupled boson-fermion mixture. If 
gBF were zero, both the bosonic and fermionic conden- 
sates would have been stable. An attractive or negative 
gBF can make one or both systems to become highly at- 
tractive and lead to collapse. Similar effect was found 
before in the study of two coupled bosonic condensates 
|l6j. In that study it was found that even if the interac- 
tion among atoms of the same species is repulsive an at- 
tractive interaction among atoms of two different species 
can cause one or both components of the coupled system 
to collapse. Similarly, in a coupled atomic-molecular con- 
densate a strong attraction between atoms and molecules 
can cause one or both of atomic and molecular conden- 
sates to collapse |l7j . 

To study the phenomenon of collapse we have to add 
the proper mechanism for atom loss in Eqs. (|2.8|l and 
(|2.9|) . The attractive interaction between a bosonic atom 
B and a fermionic atom F leads to the following three- 
body recombination process to form a molecule (BF) 



composed of a bosonic and fermionic atom, which is re- 
sponsible for the loss of atoms p| 



B + B + F -> (BF)+B. 



(2.10) 



As the fermions are in a spin-polarized state, due to Pauli 
principle two-fcrmion molecules cannot be formed. Con- 
sequently, femionic atoms could only be lost in the pres- 
ence of bosons and the loss rate scales quadratically with 
bosonic density and is independent of fermion number 
Np [3 • I n addition to reaction H2.10fl . there is also the 
possibility of the formation of a two-boson molecule (BB) 
via the reaction 



B + B + B 



(BB)+B. 



(2.11) 



However, for a low-density repulsive bosonic condensate 
as in the experiment on boson-fermion mixture |5(, the 
rate of loss of bosons due to reaction l|2.11|l will pre- 
sumably be small j^]. Moreover, reaction (|2.11l) has no 
influence on the collapse of the fermionic condensate. As 
the primary motivation of this paper is to investigate the 
collapse of the fermionic condensate, we neglect reaction 
(|2.11|) in our study. 

The recombination reaction l|2.1U|l adds an ex- 
tra term in Lagrangian density (|2.6I) of the form 
— ihK^ p\ 2 \^ b\ 4 proportional to the product of 
fermion density and the square of boson density, where 
K3 is the loss rate. In the presence of only the bosonic 
atoms the contribution to the Lagrangian density for the 
recombination reaction (|2.11|) is the usual three-body re- 
combination term — ihKy,\^\ 6 /2 0,0] proportional to 
the cube of boson density. After proper symmetrization, 
with the use of the loss term -i^ 3 |\I/ F | 2 |* B | 4 /2 in Eq. 
(E3, Eqs. (E2J) and (£3 become 
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In the spherically symmetric state of zero angular mo- 
mentum the wave function has the form ^(r; t) — ip(r; t). 
Of the experimental condensation of boson-fermion mix- 
tures the two pioneering possibilities were the mixture of 
6 Li (fermion) and 23 Na (boson) |3L and the mixture of 
40 K (fermion) and 87 Rb (boson) [3. The actual exper- 
iment on the collapse of fermions |5( was performed in 
the 40 K- 87 Rb system with an axially symmetric trap. In 
this paper we shall consider the effect of the variation of 
the boson-fermion scattering length asF on the collapse. 
For that purpose one requires not only an accurate value 
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of -K3 but a variation of K3 with clef- As the dynam- 
ics of collapse is sensitive to this unknown variation of 
-K3, we are not in a position to produce the (quantita- 
tive) dynamics of collapse in the realistic (experimental) 
situation of Ref. |5|. Instead, in this paper we consider 
a spherically-symmetric model to see if this model can 
describe the essentials of the observed collapse dynamics 
[j| . We also present results on central probability density 
of the condensate to confirm the collapse. 

As we shall not be interested in a particular boson- 
fermion system in this paper, but will be concerned with 
the collapse of the fermionic condensate in general we 
take in the rest of this paper mj = 3m p = m(Rb), 
whence ma = 3?tif/4. In the above two experimental 
situations 0, 0| mg ~ 3mp- Also, we take Vs(r) = 
VHr) = mB0J B r 2 /2 which corresponds to a reduction 
of lo b by a factor ym^/rop as in the study by Mod- 
ugno et al. |20| . These two assumptions give a simpler 
analytical form of the final equations we derive eliminat- 
ing factors of masses and frequencies without any con- 
sequence to our qualitative study. Now transforming to 
dimensionless variables defined by x — ^/2r/l, t = tu, 
I = ^jTij (mB<jJB): and 

x 

i = B,F, we get 

d d 2 
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(2.16) 



where Nbb = Nbclbb/1, Nbf = N F a BF /l, Nfb = 
N b clbf/1, £,b = K^/ (2n 2 a B BaBFl i ^B), and £ F = 
K3/ (2ir 2 o bb 1 4 lob) ■ In the nonabsorptive case £ F = £s = 
0, the normalization of the wave-function components is 
given by J °° dx\<pi(x; r)| 2 = l,i = B,F. In the absorp- 
tive case £p 7^ and £b 7^ 0, the normalization reduces 
with time due to loss of atoms. 



III. NUMERICAL RESULT 

We solve the coupled mean-ficld-hydrodynamic equa- 
tions (|2.15|) and (|2.16|) numerically using a time-iteration 
method based on the Crank-Nicholson discretization 
scheme elaborated in Ref. 22;]. We discretize the CP 
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FIG. 1: (Color online) The evolution of the numbers of bosons 
iVfl(i) and fermions N F (t) during collapse initiated by a jump 
in boson-fermion scattering length clbf from —12.5 nm to 
-37.5 nm for K 3 = 2 x 1(T 28 cm 6 s _1 , K 3 = 1(T 27 em's -1 , 
and K3 = 10 -26 cm 6 s _1 . The three upper (blue, dotted) 
curves refer to number of bosons and three lower (red, solid) 
curves refer to number of fermions. The curves are labeled by 
their respective K3 values. 



equation using time step 0.00025 and space step 0.05 
spanning x from to 25. This domain of space was suf- 
ficient to encompass the whole condensate wave function 
during collapse and expansion. 

First we solve Eqs. I|2.15|l and I2.16f> with £ B = £ F = 0. 
This will allow us to find a stable bound state of boson- 
fcrmion mixture. After some experimentation we take 
in our calculation clbb — 5 nm, clbf = —12.5 nm, 
lob = 2tt x 100 Hz, N B = 4800, and N F = 1200, 
so that I w 1 fim, o.bb/1 — 0.005, o.bf/1 = —0.0125, 
Nbb = 24, Nbf = -15, and Nfb = -60. These values 
of parameters are similar to those employed in experi- 
ments. The dimensionless unit of time corresponds to 
tug 1 ps 1.6 ms, and dimensionless unit of length corre- 
sponds to 1/V2 w 0.7 /im. 

Now we consider the collapse of fermions initiated by 
a sudden jump in the boson-fermion scattering length 
from asF = —12.5 nm to —37.5 nm which can be impli- 
mented near a boson-fermion Feshbach resonance. Such 
boson-fermion Feshbach resonances between 23 Na (bo- 
son) and 6 Li (fermion) atoms 1231 and between 87 Rb (bo- 
son) and 40 K (fermion) atoms [24| have been experimen- 
tally observed. These resonances should enable exper- 
imental control of the interspecies interactions [24| and 
hence can be used to increase the attractive force between 
bosons and fermions which in turn increases the attrac- 
tive nonlinearities A\/2Nbf and \\[2Mfb hi Eqs. (|2.15() 
and H2.16(l . If these attractive nonlinear terms are larger 
than the repulsive nonlinearities in these equations it is 
possible to have collapse of fermions or bosons or both. 
In order to have collapse, the effective nonlinearities in 
these equations should be attractive. 

Due to the imaginary terms in Eqs. (|2.15|l and (|2.16(l 
the numbers of bosons and fermions decay with time. 
When the net nonlinear attraction in these equations is 
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FIG. 2: (Color online) The initial (t = 0) and final (t = 
50 ms) fermion and boson wave functions in arbitrary units 
before and after collapse initiated by a jump in boson-fermion 
scattering length clbf from —12.5 nm to —37.5 nm. 



small there is smooth and steady decay of number of 
atoms. However, when the net nonlinear attraction is 
gradually increased the steady decay of number of atoms 
develops into a violent decay called collapse. When this 
happens the condensate loses a significant fraction of 
atoms in a small interval of time (milliseconds) after 
which a remnant condensate with a reasonably constant 
number of atoms is formed. Also, during and immedi- 
ately after collapse, the wave function of the condensate 
becomes very unsmooth and spiky in nature as opposed 
to a reasonably smooth wave function in the case of a 
steady decay. 

We study the the evolution of the boson and fermion 
numbers in the condensate from time t — to t = 50 
ms after the sudden jump in the scattering length from 
asF = —12.5 nm to —37.5 nm at time t = 0. The evolu- 
tion of boson and fermion numbers depends on the value 
of the loss rate K 3 . We study the sensitivity of the re- 
sult on K3 by performing the calculation for different loss 
rates. In Fig. 1 we plot the evolution of the boson and 
fermion numbers for different loss rates: A3 = 2 x 10 ~ 28 
cm 6 s _1 , 10~ 27 cm 6 s _1 , and 10~ 26 cm 6 s~ 1 . With the in- 
crease of A3 the rate of decay increases, although the 
results for different A 3 are qualitatively similar. We see 
in Fig. 1 that in all cases both the number of bosons 
and fermions decay rapidly and attain an approximately 
constant (remnant) number in less than 50 ms. The 
panorama is similar to the collapse in attractive bosonic 
condensate studied experimentally by Donley et al. |ll| , 
where also a cold remnant bosonic condensate is formed 
at the end of the collapse. In the recombination pro- 
cess <|2.10ll two bosonic and one fermionic atoms are lost. 
In Fig. 1 we find that during the collapse about 2000 
bosonic atoms and 1000 fermionic atoms are lost. 

Experimentally, Modugno et al. measured the follow- 
ing loss rate A 3 = 2(1) x 10~ 27 cmS" 1 for the 40 K- 87 Rb 
system |5( . In the remainder of this paper we shall use the 
value A3 = 10~ 27 cm 6 s _1 for all densities and all values 
of scattering lengths. This value of A3 gives a loss rate 
during collapse compatible with experiment [f|, where a 
significant fraction of fermions are lost in an interval of 
time A < 50 ms. We further examined the wave function 
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FIG. 3: Evolution of the (a) central boson probability density 
(Nb\i/jb(0, t)\ 2 ) and (b) central fermion probability density 
(Nf\iPf(0, t)\ 2 ) in arbitrary units during the collapse exhib- 
ited in Figs. 1 and 2. 



of the bosons and fermions to determine whether the de- 
cay in Fig. 1 really corresponds to a collapse and not to 
a slow evaporation due to the imaginary terms in Eqs. 
(l2~l"5l) and lf2~TBj) . 

In Fig. 2 we plot the profiles of the bosonic and 
fermionic wave functions in arbitrary units at times t = 
and t — 50 ms. A close look at Fig. 2 reveals that before 
collapse at t = the bosonic and fermionic wave func- 
tions are smooth and the fermionic wave function extends 
over a larger distance than the bosonic wave function. 
This is due to the large repulsion between the fermions 
in the spin polarized state due to the Pauli exclusion 
principle. Hence the bosonic condensate lies well inside 
the fermionic condensate. This was noted in experiment 
@ , as well as in previous theoretical studies El H3 • 
The wave functions after collapse have an entirely differ- 
ent profile. As expected the wave functions are highly 
peaked in the central (r = 0) region and develop spikes. 
However, they extend over a large distance too. The final 
spiky wave function indicates the collapse in contrast to 
a smooth final wave function corresponding to a steady 
loss of atoms. The collapse is a quick process lasting at 
most a few tens of milliseconds when a significant fraction 
of atoms are lost. For example, in Fig. 1 for A3 = 10~ 27 
cm 6 s -1 , the collapse lasts for the first 25 ms when most 
of the atoms are lost. After this interval the rate of loss 
of atoms is reduced and remnant bosonic and fermionic 
condensates with a roughly constant number of atoms are 
formed. The wave function after the collapse of a purely 
bosonic condensate also exhibits a similar behavior [l4j . 

To confirm further the collapse in Figs. 1 and 2 for 



A'3 = 10 27 cm 6 s 1 , we plot in Fig. 3 the evolu- 
tion of the central probability density of (a) the bosonic 
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FIG. 4: (Color online) (a) The evolution of the numbers of 
bosons NB(t) and fermions Nf(£) during collapse initiated 
by a jump in the boson number from 4800 to 14400 for K$ — 
10" 27 cm 6 s~\ (b) The initial (t = 0) and final (t = 50 ms) 
fermion and boson wave functions in arbitrary units before 
and after the collapse in (a). 



(N B \4>B(0,t)\ 2 ) and (b) fermionic {N F \ip F (0, t)\ 2 ) com- 
ponents in arbitrary units during collapse. We note a 
very strong fluctuation of the central density reminiscent 
of collapse both in the bosonic and fermionic components. 
Similar variations are common to the collapse of a purely 
bosonic condensate Such strong fluctuation of the 

central density could not be due to a weak evaporation 
of the condensate due to recombination l|2.1()Jl . 

Next we consider fermionic collapse initiated by an in- 
crease of the number of bosonic atoms. In Fig. 4 (a) we 
plot the evolution of the boson and fermion numbers in 
the condensate from time t = to t = 50 ms after the 
jump in the boson number from 4800 to 14400 at time 
t = 0. In this case there is a decay in boson and fermion 
numbers. The decay of boson numbers (~ 2000) is dou- 
ble the decay of fermion numbers (~ 1000). To determine 
if this case corresponds to collapse, in Fig. 4 (b) we plot 
the profiles of the bosonic and fermionic wave functions 
in arbitrary units at times t = and t = 50 ms. At 
t = 50 ms the fermionic wave function is sharply peaked 
at the center and spiky in nature, whereas the reasonable 
smooth bosonic wave function does not exhibit any cen- 
tral peaking. The slowly varying final bosonic component 
corresponds to a weak loss of atoms and not to collapse. 
The collapse of the fermionic component deserves further 
examination. In Fig. 5 we plot the evolution of the cen- 




FIG. 5: Evolution of the central fermion probability density 
(Nf\iPf(0, t)\ 2 ) in arbitrary units during the collapse exhib- 
ited in Figs. 4. 



tral probability density of the fermionic (N F \ipg(0,i)\ 2 ) 
component in arbitrary units during collapse exhibited 
in Figs. 4. The central density exhibits rapid oscillation 
indicating collapse. However, the fluctuation in Fig. 5 is 
less than those noted in Figs. 3 indicating a less violent 
collapse in this case. This is quite reasonable as in Fig. 
2 both components undergo collapse. 

The initial nonlinearities in Eqs. (|2.15|l and (|2.16|) are 
Nbb = 24, Nbf = —15, Nfb = —60, whereas the final 
nonlinearities in Figs. 1 are Nbb — 24, Nbf = —45, 
Nfb — —180. The final attractive (negative) nonlin- 
earities are so strong that the effective nonlinearities in 
both Eqs. I|2.15|l and l|2.16[) have become attractive and 
large in nature. This is responsible for the collapse ob- 
served in both components in Figs. 1 and 2. In the 
situation of Fig. 4 the final nonlinearities are Nbb = 72, 
Nbf = —15, Nfb = —180. We find that, as a result, the 
repulsive nonlinearity in the bosonic equation (|2.15|) has 
increased and hence the bosonic condensate does not col- 
lapse; whereas the attractive nonlinearity in the fermionic 
equation i|2.16() has increased by a factor of 3 and hence 
this component undergoes collapse. However, in Fig. 4 
(a) there is a steady loss of atoms in the repulsive bosonic 
condensate during the collapse of the fermionic conden- 
sate. 

The present study indicates that under different situa- 
tions the fermionic component in a trapped condensed 
boson-fermion mixture can undergo a collapse due to 
an attractive boson-fermion interaction. In the present 
study we consider two situations. Of these, the fermionic 
collapse associated with the increase in the number of 
bosons has been exploited in the experiment of Mod- 
ugno et al. p|- In the present model simulation, we 
introduced a sudden jump in the number of atoms which 
is not possible in the laboratory. In the experiment of 
Ref. |3 a steady increase in boson numbers has been 
used to initiate the fermionic collapse. However, the sit- 
uation of collapse exploiting a Feshbach resonance in the 
boson-fermion system is more exciting and controlled ex- 
periment could be done by jumping the boson-fermion 
scattering length clef as in the collapse experiment with 
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bosonic condensat e [Til . This seems possible in 6 Li- 23 Na 
[23| and 40 K- 87 Rb |24| systems using the recently discov- 
ered Feshbach resonances in these systems. The increase 
in the number of bosonic atoms in the boson- fermion con- 
densate to initiate the collapse is a slow stochastic process 
as opposed to sudden controlled jump in the scattering 
length leading to a much violent collapse. In this con- 
nection we recall that the collapse of a purely bosonic 
7 Li condensate in Ref. 10 by a steady increase of the 
atom number was a stochas tic p rocess, whereas the col- 
lapse of 85 Rb atoms in Ref. [Tl[ by a jump in scattering 
length was a more violent, nevertheless more exciting, 
process. However, the collapse initiated by a jump in 
the scattering length clbf will require a careful study of 
the boson-fermion system in the search of an appropriate 
Feshbach resonance [23|, |2J] . 

An examination of Eqs. I|2.15[l and l|2.16[l with £p = 
£b = reveals that for a fixed Np with the increase of 
the boson number Nb the attractive nonlincarity Mfb 
can overcome the repulsive Fermi pressure and lead to a 
collapse of the fermionic condensate. Hence, for a fixed 
Nb and fixed values of scattering lengths, Np has to 
be larger than a critical value in order to have a sta- 
ble condensate. The fermionic collapse is a function of 
the boson-fermion scattering length cibf and the boson 
number Nb ■ These aspects of collapse have been studied 
via equilibrium equations (|2.3[) and l|2.4|l by Modugno 
et al. 20] . Finally, we comment that a large number 
of fermions Np in Eqs. (|2.15|l and (|2.16|) will lead to a 
strong attraction in the bosonic equation l|2.15|l via the 
term Mbf which may initiate a collapse of the bosonic 
atoms. An increase of Nb stabilizes the bosonic con- 
densate but may initiate collapse in the fermionic one, 
whereas an increase in Np stabilizes the fermionic con- 
densate but may start collapse in the bosonic one. Fi- 
nally, an increase in |obf| may initiate collapse in both 
components. A detailed study of these features in the 
actual experimental situation would be a welcome future 
work. 



IV. SUMMARY 

We have suggested a coupled set of time-dependent 
mean-field-hydrodynamic equations for a trapped boson- 
fermion condensate using the (time-independent) energy 
functional of Ref. 9] successfully used to study the 
equilibrium states of the trapped boson-fermion conden- 
sate. One could identify collapse from the time inde- 
pendent formulation, where equilibrium stationary so- 
lutions for the boson-fermion mixture cease to exist. 
The present time-dependent generalization permits us 
to study nonequilibrium dynamics of the coupled boson- 
fermion condensate. We use the time-dependent nonlin- 
ear model to study the collapse dynamics for the boson- 
fermion mixture. There arc two possibilities for the col- 
lapse of the fermionic component. From Eq. (|2.1ti[) we see 
that this happens when the attractive nonlincarity Mf b 
becomes stronger either via an increase in boson number 
Nb or via an increase in the strength of boson-fermion 
interaction |obf|- We considered both possibilities in the 
present numerical study. The increase of |obf| leads to a 
collapse of both components, whereas the increase of Nb 
leads to a collapse of the fermionic component alone. The 
collapse is more violent in the first case. The present sim- 
ulation was done in a spherically symmetric model. The 
collapse dynamics is strongly dependent on the loss rate 
K3 . Once one has a good knowledge of the variation of 
K3 with |aBiH, an d high-quality experiment data will be 
available, a more realistic calculation will be worth doing 
in the future. 



Acknowledgments 

The work is supported in part by the CNPq of Brazil. 



[1] B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 
[2] K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. 

Granade, and J. E. Thomas, Science 298, 2179 (2002). 
[3] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. 

Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. Lett. 

87, 080403 (2001); A. G. Truscott, K. E. Strecker, W. I. 

Mc Alexander, G. B. Partridge, and R. G. Hulet, Science 

291, 2570 (2001). 
[4] Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, M. 

W. Zwierlein, A. Gorlitz, and W. Ketterle, Phys. Rev. 

Lett. 88, 160401 (2002). 
[5] G. Modugno, G. Roati, F. Riboli, F. Ferlaino, R. J. 

Brecha, and M. Inguscio, Science 297, 2240 (2002). 
[6] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 

Rev. Lett. 89, 150403 (2002). 
[7] K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. 

Rev. Lett. 91, 080406 (2003); Z. Hadzibabic, S. Gupta, 



C. A. Stan, C. H. Schunck, M. W. Zwierlein, K. Dieck- 
mann, and W. Ketterle, ibid. 91, 160401 (2003). 

[8] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[9] P. Capuzzi, A. Minguzzi, and M. P. Tosi, Phys. Rev. A 
69, 053615 (2004); 68, 033605 (2003); 67, 053605 (2003). 
[10] J. M. Gerton, D. Strekalov, I. Prodan, and R. G. Hulet, 

Nature 408, 692 (2001). 
[11] E. A. Donley, N. R. Claussen, S. L. Cornish, J. L. 
Roberts, E. A. Cornell, and C. E. Wieman, Nature 412, 
295 (2001) 

[12] S. Inouye, M. R. Andrews, J. Stenger, H. J. Miesner, 

D. M. Stamper-Kurn, and W. Ketterle, Nature 392, 151 
(1998). 

[13] S. K. Adhikari, Phys. Rev. A 66, 013611 (2002); 66, 
043601 (2002). 

[14] H. Saito and M. Ueda, Phys. Rev. A 65, 033624 (2002); 



8 



L. Santos and G. V. Shlyapnikov, ibid. 66, 011602 (2002); 
C. M. Savage, N. P. Robins, and J. J. Hope, ibid. 67, 
014304 (2003); R. A. Duine and H. T. C. Stoof, ibid. 68, 
013602 (2003); E. A. Calzetta and B. L. Hu, ibid. 68, 
043625 (2003); W. Bao, D. Jaksch, and P. A. Markowich, 
J. Phys. B 37, 329 (2004); S. K. Adhikari, ibid. 37, 1185 
(2004). 

[15] H. Heiselberg, C. J. Pethick, H. Smith, and L. Viverit, 
Phys. Rev. Lett. 85, 2418 (2000); M. J. Bijlsma, B. A. 
Heringa, and H. T. C. Stoof, Phys. Rev. A 61, 053601 
(2000). 

[16] S. K. Adhikari, Phys. Rev. A 63, 043611 (2001); Phys. 

Lett. A 281, 265 (2001). 
[17] S. K. Adhikari, J. Phys. B 34, 4231 (2001). 
[18] V. S. Shchesnovich, A. M. Kamchatnov, and R. A. 

Kraenkel, Phys. Rev. A 69 033601 (2004). 
[19] R. Roth, Phys. Rev. A 66, 013614 (2002); R. Roth and 



H. Feldmeier, ibid. 65, 021603 (2002); T. Miyakawa, T. 

Suzuki, and H. Yabu, ibid. 64, 033611 (2001); X.-J. Liu, 

M. Modugno, and H. Hu, ibid. 68, 053605 (2003); X.-J. 

Liu and H. Hu, ibid. 67, 023613 (2003). 
[20] M. Modugno, F. Ferlaino, F. Riboli, G. Roati, G. Mod- 
ugno, and M. Inguscio, Phys. Rev. A 68, 043626 (2003). 
[21] Yu. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 

Rev. Lett. 79, 2604 (1997). 
[22] S. K. Adhikari and P. Muruganandam, J. Phys. B 35, 

2831 (2002); P. Muruganandam and S. K. Adhikari, ibid. 

36, 2501 (2003). 
[23] C. A. Stan, M. W. Zwierlein, C. H. Schunck, S. M. F. 

Raupach, and W. Ketterle, cond-mat/0406129 
[24] S. Inouye, J. Goldw in, M. L. Olsen, C . Ticknor, J. L. 

Bohn, and D. S. Jin, \cond-m at/0406208 



